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in electrostatically driven granular layers. 

PACS numbers: 05.40.-a, 05.20.Dd, 5.45.-a 



I. INTRODUCTION 

A multitude of growth phenomena in physical pro- 
cesses are driven by exchange of particles between clus- 
ters. Examples include droplet growth via evaporation 
and re-condensation y| , island growth in deposition pro- 
cesses , and phase ordering [ESQ- Exchange pro- 
cesses have been also used to model social and economi- 
cal systems including segregation of heterogeneous popu- 
lations 0, the distribution of wealth in a society Q, and 
growth of urban populations @, . 

In exchange processes, clusters are composed of 'atoms' 
(monomers). Monomers detach from one cluster and re- 
attach to another cluster. We shall consider the detach- 
ment controlled limit where the time scale for transport 
between clusters is much faster then the time scale for 
detachment. Exchange processes incorporate both re- 
versible and irreversible features. Clusters may grow or 
shrink, yet when a monomer attaches to another clus- 
ter, its respective cluster disappears. This irreversible 
step provides the mechanism for cluster growth. There- 
fore, exchange-driven processes are fundamentally differ- 
ent from irreversible growth processes, particularly ag- 
gregation in mm. 

Such mass transfer processes are governed by an ex- 
change kernel K(i,j) that represents the rate of transfer 
of monomers from a cluster of size i to a cluster of size 
j. Generally, the rate of monomer exchange between two 
clusters depends on their sizes. Moreover, we consider the 
case where there is no preferable direction for exchanges, 
i.e., symmetric exchange kernels, K(i,j) — K(j,i). This 
is unlike migration processes where the exchange is pref- 
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erential ("big gets bigger" or "rich gets richer"). Migra- 
tion underlies certain physical processes (e.g., coarsening 
with conserved order parameter 0,3) as well as social 
and economical processes 0, IS Et ■ 

We investigate homogeneous exchange kernels, 
K(ai,aj) — a 2X K(i,j). In particular, we consider the 
product kernel K(i,j) = and its generalization 

K(i,j) = + iPj" with v + fi = 2A and v > fi. 
We obtain a complete description of the problem in 
the asymptotic scaling regime. The overall range of 
possible behaviors and the emergence of self-similar size 
distributions are as in aggregation and migration pro- 
cesses. However, there are quantitative and qualitative 
differences. Unlike aggregation, the gelation transition 
is complete, and unlike migration, the size distributions 
are extended rather than compact. 

The behavior falls into three categories. I) Growth: 
When v < 2 and A < 3/2, clusters grow indefinitely. 
The typical cluster size grows algebraically with time, 
k ~ t 1 /^ 3_2A \ and the cluster size distribution is given 
by a self-similar distribution with a stretched exponential 
tail. II) Gelation: When v < 2 and A > 3/2, the entire 
mass in the system is transformed into an infinite gel 
in a finite time. The cluster mass diverges algebraically 
near the gelation point, k ~ (t c — t) 1 /( 3_2A ) and a scal- 
ing behavior similar to the one underlying the growth 
phase is found. In the borderline case v = 2 the scaling 
function has an algebraic tail with a universal exponent 
<$>(x) ~ £~ 5 . Scaling breaks down in the special point 
v = fi = 2 where the distribution is log-normal. Ill) In- 
stant gelation: When v > 2, the gelation time vanishes 
logarithmically with the system size, t c ~ [In A^] - ^ -2 ). 
In particular, for an infinite system, gelation is instanta- 
neous. 

This paper is organized as follows. In the next section, 
we define the exchange process. The governing equations 
are analyzed using scaling techniques and exact solutions 
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for the moments. We first analyze the product kernel 
(Sec. Ill) and then, the generalized kernel (Sec. IV). The 
gelation time in finite systems is investigated in Sec. V us- 
ing heuristic arguments and numerical simulations. Ap- 
plications to coarsening in the Ising model with infinite 
range Kawasaki dynamics and in electrostatically driven 
granular layers are briefly discussed in Sec. VI, and con- 
clusions are given in Sec. VII. 



II. EXCHANGE PROCESSES 

We consider the following elementary exchange pro- 
cess. The system consists of an ensemble of clusters and 
clusters evolve via transfer of a single monomer from one 
cluster to another. Symbolically, 



{i,j) K M ] (iil.jTl) 



(1) 



with i and j the number of particles in each cluster and 
K(i,j) the exchange kernel. In an exchange event, a 
cluster is equally likely to gain or to lose a particle. Since 
the exchange process is unbiased, the matrix of transition 
rates is symmetric: K(i,j) = K(j,i). 

Let A k (t) be the density of clusters containing k 
monomers at time t. It evolves according to the following 
rate equation 



dA 



—tt = ^2 AiA-jK(i,j) [Sk,i+i + §k,i-i ~ 2<5fc,i 



(2) 



This equation assumes perfect mixing, or equivalcntly, 
absence of spatial correlations. We restrict our attention 
to monodisperse initial conditions, A k (0) = 8 k> \. The ex- 
change process has a single conservation law. As reflected 
by the evolution equations, the total mass is conserved, 
Mi = 1 with M a = ^2 n k a Ak(t) the moments of the 
size distribution. It is natural to consider homogeneous 
kernels, K(ai,aj) — a 2X K(i,j) with 2A the homogene- 
ity degree and we present results for the product kernel 
K(i,j) — (ij) x and the generalized homogeneous kernel 
K(i,j) = + i^f with \i + v = 2A. 



III. THE PRODUCT KERNEL 

For the product kernel, K(i,j) = (ij) x , the rate equa- 
tions read 



dA k 
dt 



= M\ [(k + l) x A k+1 + (k- l) A A fe _! - 2k x A k ] 



with the boundary condition Aq = 0. These evolution 
equations demonstrate the diffusive character of the ex- 
change process. Absorbing the factor M\ into the time 
variable 



t= dt' M x (t') 
Jo 



(3) 



we recast the governing equations into 
dA k 



dr 



(k + l) x A k+1 + (k - l) A ^_i - 2k x A k 



(4) 



Alternatively, one can study integer moments of the size 
distribution. The total density obeys -§^Mq = —A\, the 
total mass is conserved -j^Mi — 0, and higher integer 
moments satisfy the following hierarchy of equations 
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dr 



i=i 
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Only for integer values of the homogeneity index is this 
hierarchy closed. We employ different approaches for dif- 
ferent A's. For A < 2, we perform a scaling analysis of 
the rate equations and for A > 2, we analyze the moment 
equations. This general analysis is augmented by exact 
solutions for the integer values A = 0, 1, and 2. 



A. Scaling (A < 2) 

When A < 2, dimensional analysis of Eq. (0} shows 
that the typical cluster size grows as 



with 



1 



2 - A 



(G) 



Using £l — M\ ~ T a< - X ^ , the growth of the typical scale 
is expressed in terms of the physical time 



' t 13 A < 3/2, 

exp(const. xt) A = 3/2, 



(t - ty 



3/2 < A < 2. 



(7) 



The dynamical exponent is (3 — (3 — 2A) _1 . As long as 
A < 3/2, clusters grow indefinitely and the characteristic 
size grows algebraically with time. For A > 3/2, a gela- 
tion transition occurs, i.e., the system develops a giant 
cluster in a finite time t c . 

We seek a scaling solution of the rate equations 



A k (r) 



-2a 



$ (k- 



(8) 



Mass conservation dictates the normalization J\ = 1 
where J a = j dx x a $(x) is the a th moment of the scaling 
distribution. Technically, the scaling function describes 
the behavior in the limits k — > oo, r — > oo with the vari- 
able x = kr~ a fixed. Thus, we consider the continuum 
limit of the rate equation -^A(k,T) — Jp- [fc A A(fc,r)]. 
The scaling function satisfies the second order linear dif- 
ferential equation 



(2 - A) ^ [x x $(x)] +z^Uh) 



dx 



2$(x) = 0. (9) 



Multiplying this equation by x, employing the identities 
x 2 & + 2x$ = (x 2 $)', x*" = (x^)" - 2*', and integrat- 
ing once yields (2 - A) [(x x+1 $)' - 2x A $] + x 2 <&(x) = 0. 
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Integrating a second time gives the scaling function: 

„2-A 1 

1-A . 



= CV~ A exp 



(2-\y 



(10) 



B 2 = 2B\. Solving recursively yields B k = fc-Bi = kA\ 
or A* k = k 1 ~ x A* 1 . Mass conservation, ^ k kA* k = 1, im- 
plies A\ = 0, and thence A* k = for all fc, i.e., complete 
gelation. 



with C = (2 - A)- 2 /( 2 ~ A Vr ^^M found from the con- 
dition Ji = 1 . The nature of the scaling function differs 
from that found for migration where K(l,m) — for 
I < m 8]: Exchange is characterized by extended dis- 
tributions, while migration is characterized by compact 
distributions. 

There are two physically relevant cases for which the 
rate equations can be solved exactly. When the exchange 
kernel is independent of the cluster size (A = 0), the rate 
equation is -^A k = A k+ i + A k ^\ — 2A k and the cluster 
size distribution is Q 



A k 



-It 



[/ fe _ 1 (2r)-7 fc+1 (2r)] 



(11) 



where /„ is the modified Bessel function of order n |13| . 
In agreement with the general scaling analysis, the typ- 
ical scale grows diffusively, k ~ r 1 / 2 , and the scaling 
function is given by $(2;) = (47r) _1 / 2 x exp(— ir 2 /4). 

For the pure product kernel (A = 1), the rate equa- 
tions read j- t A k = (k + l)A k+ i + (k - l)A fc _i - 2kA k 
(in this case t = r). Substituting the mass-conserving 
ansatz A k = (1 — u) 2 u k ~ 1 reduces the infinite set of 
rate equations into a single ordinary differential equation 



dt 

The size distribution in this case is 



u 



(1 — u) subject the initial condition u(0) = 0. 



(1 + t) k + 1 



(12) 



The typical cluster size grows ballistically, k ~ t, and the 
scaling function is purely exponential, $(2;) = e~ x , again 
in agreement with the above scaling results. 

When 3/2<A<2, an infinite cluster is formed at 
some finite time t c , termed the gelation time. The gela- 
tion time depends on the initial condition and its de- 
termination requires the full time dependent behavior. 
Even without knowing the gelation time exactly, one can 
describe the behavior in the pregel stage since the size 
distribution still admits the scaling form (|1U|) . Thus, for 
all A < 2 we have 



A k {r) -Ck 1 



exp 



(2 -A) 2 



(13) 



implying that A k — > in the limit r — > 00 (t — > t c ). 
In other words, the gelation is complete, A k (t) — 0, for 
t > t c . This behavior differs from aggregation processes 
where at the gelation time, all mass in the system is 
contained in finite size clusters 0, 0, 0] . 

Complete gelation can be alternatively shown as fol- 
lows. Let us assume that the cluster size distribution 



approaches a constant A k 
Eqs. (0J, the quantities B k = 
Laplace equation B k+ \ + B k _\ 



At > as r — » 00. From 



k x A k satisfy the discrete 



x k 
2B k 



for k > 1 and 



B. Multiscaling (A = 2) 

In this special case, the moment equations (JSJ) are 
closed for n > 2. For example 



dM 2 

~dV 

dM 3 

~~dT~ 

dM 4 

~dV 



= 2M 2 , 
= 6M 3 , 

= 12M 4 + 2M 2 . 



(14) 



The solutions to these equations are combinations of ex- 



ponentials: M 2 



„2t 



.fir 



etc. The physical time t - ., 
t = J Q T dsM 2 ' 1 (s) so 



|[1 



-2t 



] is found from 



M 2 = (1-2*)- 1 , 
M 3 = (l-2t)-\ 

M 4 



(15) 



\(l-2t)-"-\{l-2t)-\ 



Therefore, the gelation time is t c = 1/2. Asymptotically, 
the first term in JSJ) dominates: ^M n ~ n(n — l)M n im- 
plying M n ~ exp[n(n — l)r] for n > 1. Close to the 
gelation time (t — * t c ), the moments diverge according to 



Af„(i)~(t c -tr n( "- i)/2 . 



(16) 



In this special case, moments exhibit multiscaling 
asymptotic behavior, namely, the normalized moments 
Mn /Mi diverge. 

To determine the asymptotic form of the size distribu- 
tion we treat k as a continuous variable. For A = 2, 
Eq. becomes -§=A k = ■^ I [k 2 A k \. This equation is 
equidimensional in k [Lj thereby suggesting use of the 
variable £ = In k instead of k. Making the transformation 
from A k (t) to A(£, r) defined via A k dk = A(£)d£, we re- 
cast above equation for A k (t) into the following constant 
coefficients diffusion-convection equation 



0_ 

7h 



d_ 



9_ 2 

01; 



(17) 



With the initial conditions A(£,0) = S(£), the solution 
reads A(£,r) = {Attt)- 1 / 2 exp[-(£ + r) 2 /(4T)]. The orig- 
inal distribution A k — k~ 1 A(£) is log-normal 



Mr) 



(47rr)- 1 / 2 e- r / 4 fc- 3 / 2 exp 



(ln/c) 5 



4r 



(18) 



Again, the distribution vanishes at the transition point, 
i.e., the gelation transition is complete. Moreover, the 
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mass distribution is algebraic, Ak(t) ~ Mo(t) k 3//2 for 



sufficiently small masses, k <C yln j-zft- The total den- 
sity vanishes quite slowly near the transition point 

fo^-fl-aWln^) V2 • (19) 



Mo 



We note that the density follows a different law than the 
one characterizing higher than first moments (|f 6f) . 

The size distribution does not follow a scaling behav- 
ior asymptotically and the log-normal distribution is re- 
sponsible for the multiscaling behavior (| 1 fc>|) of the mo- 
ments. This differs from aggregation processes where the 
moments diverge as M n (t) oc (t c — i) -Q " 0] with the 
exponent a n linear in n. 



C. Instant Gelation (A > 2) 

Gelation is now instantaneous and complete, that is 
A k (t) = for all k when t > 0. To prove this assertion 
we assume the opposite and arrive at a contradiction. 
Our analysis follows van Dongen's approach developed 
in the context of aggregation processes . 

The moments M n with integer n > 2 evolve according 
to (Jsj. The first term in the summation yields a lower 
bound for their growth rate 



dMn 
dr 



> n(n - 1)M„_ 2 +A > n(n - 1)(M„) 



l+A 



(20) 



with A = ^Ej • The second inequality follows from the 
Jensen's inequality as shown below. Consider the auxil- 
iary functions M n , evolving according to 



dM n 
dr 



n{n-l)(M n ) 



l+A 



(21) 



Solving this equation subject to the initial condition 
M n {0) = 1 yields M n = [1 - n(X - 2)r]" 1 / A . Therefore, 
M n — > oo as t — > r n = n _1 (A — 2) -1 . Since M n > .M n , 
the moment M n diverges at least at r„. The series of 
times T n set an upper bound for the gelation time r c 
since all moments should be finite for t < r c . As t„ — » 
when n — ► oo, we conclude that r c = and thence, the 
gelation time vanishes t c = 0. 

The inequality M„_ 2 +a > (M„) 1+A with A = £=| is 
derived as follows. Let the parameters pj > satisfy 
^2jPj — 1 and let <fr(x) be a convex function. A convex 
function satisfies the Jensen inequality 



OO / OC \ 

^ /, ,<!-:.,■ , : 'i>( X /',••'•,•)• 



(22) 



First, we substitute the coefficients pj — jAj (from 
mass conservation y\ jAj = 1) and the convex function 



$(z) 



_ ^l+A 



(A > for A > 2) into the Jensen in- 



equality. Then, choosing Xj = j n 
Y..i h -\, = M n and >»:.rj 
above inequality. 



and using ^2pjZ 



M, 



n-2+\ 



yields the 



IV. GENERALIZED KERNELS 

The rates K(i,j) underlying exchange processes are 
typically homogeneous functions of i and j (at least 
for large i and j). We restrict ourselves to such ker- 
nels. Apart from the homogeneity degree 2A, homo- 
geneous kernels are characterized by an additional ex- 
ponent v defined through the asymptotic K(l,j) ~ f 
as j 3> 1. For j < j the exchange kernel scales as 



;2A 



K(l,j/i) ~ i"j u with 2A = v + fi. There- 



fore, we consider a specific generalization of the product 
kernel that exhibits these homogeneity properties 



K(i,j) = i u j> i + i> l f 



(23) 



More precisely, the asymptotics K(i,j) ~ i^j" 
for i -C j if v > (J-', since the kernel is symmetric, we 
can assume that v > [i without loss of generality. We 
expect that the homogeneity indices govern the overall 
qualitative behavior (growth, gelation, instant gelation), 
while the precise form of the kernel controls quantitative 
characteristics such as the size distribution. 

For this exchange kernel, the rate equations J2J) become 

A A. 

-± = M n [(k + l) v A k+l + (k - l)"Afc_i - 2k v A k ] 
+ M v [(k + lTA k+1 + (k- l)"Afc_i - 2k»A k ] . 

The following generalization of the modified time variable 

dt'^AU(t')M^t') (24) 

handles the two indices symmetrically. In terms of this 
time, the evolution equations are 



dA k 
dr 



= R [(k + l) v A k+1 + (k - l) v A fc _i - 2k" A k ] 
+ R- 1 [{k + iy i A k+1 + {k- lfA k ^ - 2k»A k ] 



with R = \J M^/M v . Of course, the dynamics conserve 
mass: 4-M\ = 0. Higher integer moments evolve accord- 
ing to 



dM n 
dr 



V2] 



1=1 



[RM n - 



.2l +v + R^Mn-U+n] ■ (25) 



When v < 2, the scaling analysis follows closely the 
product kernel case. The overall growth laws © and 
remain unchanged and the homogeneity degree A charac- 
terizes the scaling behavior. However, we shall see that 
the individual indices v and fi play an important role 
since they dictate the range for which this law holds. 

We seek a scaling solution of the form ijHJ. The scaling 
function $(x) satisfies 



^2 r 



dx 2 



(Ux u + Vx> M )$(x) 



x— $(x) + 2<J>(x) = (26) 
dx 
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with the constants U — a A and V — a A~ l deter- 
mined from the ratio A = \J J^/ J v . The scaling function 
reads 



$(s) = C 



Ux v + Vx^ 



exp 



dy 

o Uy" + Vy» 



(27) 



The scaling solution involves three parameters U, V, 
and C. Substituting U — a~ 1 A and V = a^ x A^ x into 
the equality A — \J J^/Jv yields a closed equation for the 
parameter A. Once A is determined, the parameters U 
and V follow, and finally, the amplitude C is found from 
the normalization J% = 1. 

We illustrate this procedure for the special case 
(v, n) = (1, 0), i.e., for the pure sum kernel K(i,j) = i+j. 
In this case, the integral on the right-hand side of Eq. 127(1 
is readily computed. Using U = ^ and V = ^ we ar- 
rive at 



$(i)=Ci(l + A 2 i)" 1 exp [-aA 2 i] 



(28) 



with a 



\A- 



We now substitute this solution 



into the right-hand side of the equality A = WJq/Ji 
and transform it into the transcendental equation 
(^) a F(a, a) + a~ 1 = 1 involving the incomplete gamma 
function (see Appendix A). The amplitude is then ex- 
plicitly evaluated to give C = aA 6 . Numerically, we find 
a ~ 2.82649, A = 0.428397, and C 0.0174713. Inter- 
estingly, there is a nontrivial algebraic correction to the 
leading exponential behavior, $(a;) ~ x a exp(— aAx) for 
large x. 

On the boundary v = 2 separating regime III from the 
two other regimes, the solution of Eq. I|27|l significantly 
simplifies. We find A = l/[2(2 - //)], U = 1/4, V = 
1/(2 — /i) 2 , and consequently, the scaling function 



$(x) = Cx 



P 2-M 



4(2 -tf 



■ (29) 

* 2 — /j ' 2 — [a n 

expressed in terms of the beta function. Remarkably, the 
scaling function l|29l) exhibits a universal large- X asymp- 
totic behavior 



The constant C = 2 [2(2 - //)]" 



[■B(, 



(30) 



In other words, the size distribution is algebraic, Ak(r) ~ 
fc" 5 r 3Q where a = (2 - A)" 1 = 2/(2 - /z). With this al- 
gebraic divergence, sufficiently small moments are char- 
acterized by ordinary scaling behavior while higher mo- 
ments exhibit multiscaling behavior: 



Mr, 



T a(n~l) 
_fm(n-l) /4 



n < 4, 
n > 4. 



(31) 



This behavior follows from the leading term in the mo- 
ment equation lj2"5|> . viz. 4-M n = n(n — l)M n R. With 



R = y/Mp/M 2 

comes £pM n = 
behavior l|31|l . 



~ At 1 and A — a/4, this equation be- 
^j—^-M n , leading to the multiscaling 




FIG. 1: The three types of behaviors: Scaling (I), Ordinary 
gelation (II), and Instant gelation (III). 



The determination of A in the general situation re- 
quires numerical evaluation, yet the form and nature of 
the size distribution is clear. For example, the minimal 
(maximal) index governs the distribution of small (large) 
clusters. Indeed, from Eq. i|27|) . the extremal behaviors 



$(x) 



x 1 -^ x < 1, 

exp[— x 2 ~ v \ x » 1. 



(32) 



Apart from the point (y, fi) — (2,2), the scaling solution 
holds for all v < 2. As in the product kernel case, growth 
occurs when A < 3/2 and gelation occurs when 3/2 < A < 
2. 

For v > 2, the scaling solution l|27|) predicts $ ~ x 1 ~ l/ . 
Such behavior is inconsistent since the moment J M di- 
verges and instead, instantaneous gelation occurs. The 
moments M n with n > 1 satisfy Eq. I|25|) and the first 
term in the summation yields a lower bound for the mo- 
ment growth 4-M n > Rn(n — l)M n _2+„. Keeping only 
this term and absorbing the factor R into the time vari- 
able, the previous proof applies. Thus, gelation is instan- 
taneous. 

To summarize, there are 3 types of behaviors, deter- 
mined by the homogeneity degrees /i and v (Fig. 1): 

I Growth. The cluster size grows indefinitely, and the 
size distribution obeys scaling. 

II Gelation. The cluster size diverges in a finite time 
and the size distribution follows a scaling solution 
near the gelation time. Gelation is complete. 

Ill Instant Gelation. The cluster size distribution van- 
ishes for alH > 0. 

The cluster size distribution exhibits a scaling behav- 
ior in regimes I and II. Scaling behavior underlies the 



6 



system everywhere except for regime III and the point 
(2,2). In the bulk of regimes I and II the size distribu- 
tion is a stretched exponential, while in the boundary 
with region III, the cluster size distribution has an alge- 
braic tail. Last, at the point (/i, v) = (2, 2) scaling breaks 
down and the distribution is log-normal. 



V. THE GELATION TIME 

Instantaneous gelation is certainly counter intuitive: 
A finite time singularity that occurs at time t = +0 ! 
Instantaneous gelation was investigated exclusively in the 
context of aggregation El El El El III Ell ■ For 
infinite systems, it is impossible to quantify the difference 
between two instant gelling. Finite systems, on the other 
hand, naturally quantify how fast a system gels. 

Consider a system consisting initially of A^ monomers. 
In a finite time tjy, all mass in the system condenses 
into a single 'runaway' cluster. How does the average 
gelation time T/v = (tjv) depend on N? When growth or 
ordinary gelation occurs, the answer follows from our pre- 
vious analysis. In the scaling regime, the growth law (J7J 
indicates that the condensation time grows algebraically 
with the system size, T/v ~ N 1 ^. In the case of ordi- 
nary gelation, the average gelation time saturates at an 
A^-independent value: T/v — > t c . The interesting case is 
instant gelation where the gelation time vanishes in the 
thermodynamic limit, T/v — > as A — > oo. 

For simplicity, we discuss the product kernel. The van- 
ishing gelation time is ultimately related to the short time 
behavior. Early on, loss terms in the rate equation Q 
are negligible and to leading order jjjAj = (j + l) A A,- + i, 
where we tacitly assumed r = t. For the initial condition 
A?'(0) = <5j,i , the leading order behavior of the density is 



(33) 



In a finite system consisting initially of A monomers, 
a j-mer first appears at time tj « jV --1 ' 3 ', 
estimated from NAjitj) — 1. For example, the first 
dimer and trimer appear at times t% = A^ 1 and tj, — 
2~( A ~ 1 '/ 2 A^~ 1 / 2 , respectively. By definition, the times 
increase monotonically, tj+\ > tj, yet the above esti- 
mates increase monotonically only for sufficiently small 
j < j*. From tj r = tj r+ i, we obtain the extremum 
(A — l) _1 lnAr using the Stirling formula. The cor- 



J* 

responding time T* = tj 
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FIG. 2: The system size dependence of the gelation time. 
Shown are the average gelation time Tjv and the normalized 
variance crjv versus the system size. The Monte Carlo simu- 
lation results correspond to an average over 10 3 independent 
realizations of the exchange process with A = 3. 



by focusing on the leading cluster that eventually grows 
to be the gel. Since it exchanges monomers back and 
forth with other clusters its growth mechanism is diffu- 
sive. For an ordinary diffusive process, 4r (fc) — 0, while 
^ (fc 2 ) = D. In our case, TJ = fc A with the typical size 
k 2 = (fc 2 ). Therefore, the typical size of the runaway 
cluster grows according to 



rffc 
~dt 



= k 



A-l 



(35) 



Integrating this rate equation from the critical size fc = 
to the system size fc = N gives the gelation time 



A 



f 



f 



■A-2 



J* 



(36) 



Since = (A — I)~ 1 lnA', the duration of the latter 
growth phase is much larger than that of the nucleation 
phase, T/v 3> T*. Therefore, the gelation time vanishes 
logarithmically 



is 



T, 



(In AO 



-(A-2) 



(37) 



T, 



A - I 
In A 



A-l 



(34) 



For later times, f > T,, the rate equations should be 
modified to account for the finiteness of the system (see 
e -g- |22l 123^1 since significant statistical fluctuations are 
induced by large runaway clusters that take over (even- 
tually only one such cluster remains) . The critical size of 
such clusters is j» ~ (A— I) _1 lnA^. As a complete analyt- 
ical solution seems out of reach, we proceed heuristically 



in the thermodynamic limit. A straightforward extension 
of the above argument to the generalized exchange kernel 
(ESI gives T N - (lnAO-^- 2 '. 

Therefore, in a finite system it may be difficult to dis- 
tinguish instantaneous gelation from the ordinary one. 
We verified the logarithmic law (|37|l numerically for 
A = 3 (Fig. 2). Probing fluctuations in the gelation 
time, we observe that the normalized variance ajj = 
(t%)/(tjs[) 2 — 1 vanishes logarithmically in the thermo- 
dynamic limit (Fig. 2). The distribution of normalized 
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gelation times becomes trivial P(i jv/Tjv) — > <5(z — 1) im- 1 
plying that the gelation time is a self-averaging quantity. 

We also examined the gelation time in two other 
growth processes, namely aggregation 

Ed, El 13 

and 

addition p3 , l25| . The above heuristic picture yields a 
similar logarithmic law albeit with a different exponent q g 

pll |. Self-averaging is observed numerically as well, and *P(x) 
we conclude that the behavior found for exchange pro- ' 
cesses is generic. 0.4 



VI. APPLICATIONS TO COARSENING 

In exchange processes, a monomer detaches from a 
cluster and subsequently re-attaches to another clus- 
ter. This elementary mechanism underlies a number of 
growth and coarsening processes. We apply our general 
theory to two coarsening processes. 

A. Infinite range Ising-Kawasaki model 

The zero temperature Ising model with infinite range 
Kawasaki dynamics evolves via exchan ges o f s p ins of op- 
posite signs belonging to domain walls |2(j l27l l2Sj . This 
model is equivalent to the Shelling's segregation model 
. In the limit of a vanishing volume fraction of one of 
the two phases, the domains of the minority phase are 
isolated and the process is essentially an exchange pro- 
cess with a product kernel K(i,j) = (ij) x . The depen- 
dence of the number of exchange candidates (i.e., spins 
in domain walls that can lower their energy by hoping 
to a different cluster) on the cluster size dictates the ho- 
mogeneity degree. For spherical clusters, only perimeter 
spins may exchange. Since the island size and the sur- 
face size grow with the radius according to k ~ R d and 
a ~ ~ fc( d ~ 1 )/ d , respectively, one has A = (d—l)/d 

and j3 = 1/(3 - 2A) = d/(d + 2). The dynamical expo- 
nent (defined through R ~ t z ) is therefore z = l/(d + 2). 
If the islands are polygons, a distinct possibility on a 
lattice, then only corner spins are active so A = 0, and 
consequently /3 = 1/3 and z = l/(3d). Both estimates 
agree in one dimension, consistent with the exact solution 
z = 1/3 0. 

B. Coarsening of thin granular layers 

In electrostatically driven granular layers, clusters nu- 
cleate around large grains [30|. Then, as charged grains 
oscillate back and forth between the two bounding plates 
due to the oscillating electric field, they may scatter of 
the plate or collide with other particle. Consequently, 
individual grains may transfer from one cluster to an- 
other. The rate of hoping into and out of a cluster is 
proportional to its area. Therefore, the homogeneity de- 
gree is unity, A = 1, implying f3 = 1 and a dynamical 
exponent of z — l/d. In two dimensions, this prediction 




X 

FIG. 3: The cumulative cluster area distribution, defined via 
^(x) = dy&(y) versus the normalized area x. The ex- 
perimental distribution represents roughly 10 3 clusters ob- 
tained from 20 different snapshots during a single realiza- 
tion of the coarsening process in which the total number 
of clusters evolved from 200 to 10. Once the average area 
is set to unity, the cumulative area distributions at differ- 
ent snapshots are the same. The theoretical distribution 
^f(x) = $(2;) = exp(— x) follows from Eq. (THfll . 



is consistent with the experimental observations z = 1/2 
[oOf. The corresponding size distribution $(ai) = e~ x 
provides a reasonable approximation for the experimental 
size distribution, obtained from a relatively small num- 
ber of clusters (see Fig. 3). The normalized variance, 
defined via a 2 = J 2 / Jf — I, is experimentally found to 
be — 0.80 ± 0.05 compared with the theoretical value 
a = l. Comparing "spatial" exchange processes where 
grains are exchanged only between neighboring clusters 
with further experimental data may help elucidate the 
relevance of spatial correlations. 

VII. CONCLUSIONS 

We have shown that kinetics of exchange processes are 
classified by the homogeneity indices of the governing 
rates. There are three possible regimes including indef- 
inite growth, gelation in a finite time, and instant gela- 
tion. Scaling behavior underlies the first two regimes. 
The size distributions are generally extended, decaying 
exponentially or algebraically for large sizes, in contrast 
with migration processes. 

We also studied the gelation time in finite systems and 
found that it decays rather slowly, following an inverse 
logarithmic law. It would be interesting to determine 
the full time dependent behavior of moments of the size 
distribution in the instant gelation regime. In the realm 
of aggregation, instantaneous gelation is relevant in as- 
trophysics, e.g., to the coalescence of planetesimals to 
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form planets or of stars to form super-massive black holes 
(see |21j and references therein). The physical relevance 
of instant gelation in exchange processes has yet to be 
demonstrated. 

Our description was on a mean-field level where all 
pairs of clusters in the system are equally likely to in- 
teract. It will be interesting to incorporate spatial fluc- 
tuations into this description. The nature of the spa- 
tial fluctuations depends on the mechanism for trans- 
porting monomers from one cluster to the other. For 
diffusive transport, one can incorporate effective fluxes 



into clusters, using the standard techniques developed 
for reaction-diffusion processes. 
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APPENDIX A: THE CASE {v,fi) = (1,0) 

Substituting Eq. ||2HJ) into A 2 = J /Ji yields 

J™dxx 2 (l + x) b - 1 e~ ax 
/ °° dxx (1 + xf- 1 e~ ax 

with b — a. Evaluation of the ratio of the integrals is 
performed as follows 



da 

= -—In 

da 



dx x (1 -h x) 



6—1 g — ax 



a—b 



--( / dx(l + x) 



b—1 g~ ax 



d 

da 



In 



da 



T(M) 



a—b 



a~ 2 + e a a~ 



r(a, a) 



